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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  MEMORANDUM  1312 


CALCULATION  OF  THE  BENDING  STRESSES  IN 
HELICOPTER  ROTOR  BLADES* 

By  P.  de  Guillenchmidt 


INTRODUCTION 


The  problem  of  determining  the  static  and  dynamic  stresses  on  the 
blades  of  rotary-wing  aircraft  in  forward  flight  has  occupied  the  atten¬ 
tion  of  many  engineers  of  every  country  for  a  long  time.  Up  to  within 
the  last  few  years,  however,  no  satisfactory  solution  has  been  found  for 
blades  having  a  distribution  of  mass  and  rigidity  varying  along  the  span. 
In  France,  Mr.  Dorand  had  for  many  years  (1932  with  the  Breguet-Dorand 
gyroplane)  used  a  graphical  method  which  gave  satisfactory  results;  how¬ 
ever,  it  was  too  long  and  called  for  precise  diagrams  on  account  of  the 
graphic  double  derivations  involved.  More  recently,  American  engineers 
have  developed  several  methods  which,  while  affording  correct  solutions 
in  the  general  case,  also  require  a  volume  of  calculation,  which  increases 
at  an  appalling  rate  when  the  number  of  points  examined  on  the  blade  are 
to  be  increased  or  higher  harmonics  for  the  external  forces  acting  on  the 
blade  are  to  be  introduced.  Accordingly,  it  has  been  necessary  to  intro¬ 
duce  more  approximate  methods,  and  it  is  these  methods  which  are  usually 
employed  in  design. 


The  purpose  of  the  present  report  is  to  describe  a  comparatively 
rapid  method  of  calculation  which  gives  a  correct  theoretical  solution 
of  the  problem  in  the  most  general  case.  This  method  is  the  result  of 
collaboration  between  the  Bureau  of  Calculation  of  the  Helicopter  Division 
of  the  National  Societies  of  Airplane  Construction  for  South  Eastern  and 
for  Central  France,  set  up  within  the  Committee  of  Rotating  Wing  Units 
of  the  French  Association  of  Aeronautical  Engineers  and  Technicians, 

(A.F. I.T.A. ) ,  on  the  instigation  of  Col.  Garry,  Chief  of  the  Section  of 
Rotating  Wing  Units  of  the  Technical  Service  Division  of  the  Air  Ministry. 

The  method  is  based  on  the  analysis  of  the  properties  of  a  vibrating 
beam,  and  its  uniqueness  lies  in  the  simple  solution  of  the  differential 
equation  which  governs  the  motions  of  the  bent  blade.  It  is  applicable, 
whatever  the  limiting  conditions  may  be  (blades  hinged,  blades  clamped, 
blades  fitted  with  dampers,  etc.  .  .  .).  It  requires,  strictly  speaking, 
the  preliminary  calculation  of  the  natural  frequencies  and  mode  shapes  of 

*  Calcul  en  Flexion  de  Pales  de  Giravions." S.N.C.A.C.  Report, 
Document  He3-0.03,  December  23,  19^. 
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the  blade  in  rotation.  This  calculation  can  be  reduced,  however,  as  will 
be  shown  later,  in  a  certain  number  of  cases,  to  the  calculation  of  the 
natural  frequency  of  the  first  model,  which  reduces  the  calculation 
required  to  some  extent. 

For  the  explanation  of  the  method,  let  us  take  the  case  of  a  hinged 
blade  in  flapping  motion  and  impose  the  usual  restrictive  asumptions, 
which  are: 

(a)  The  twisting  deflections  of  the  blade  are  negligible. 

(b)  The  blade  is  rigid  in  its  plane,  that  is  to  say,  the  drag 
deflections  are  negligible  with  respect  to  the  deflections 
normal  to  the  plane  of  the  blade. 

(c)  The  deflections  and  the  flapping  angles  are  small,  hence  their 
higher  powers  can  be  disregarded  and  we  may  assume 

cos  (3  ~  1  sin  3  ~  tan  3  ~  3 

(d)  The  bending  deflections  of  the  blade  do  not  influence  the  aero¬ 
dynamic  forces  acting  on  the  blade.  Included,  however,  is  a 
term  for  the  aerodynamic  damping  due  to  the  fact  that  the  vibra¬ 
tion  of  the  blade  produces,  for  each  element,  a  change  in  the 
relative  velocity,  and  consequently -©f  the  angle  of  attack.2 


SYMBOLS 


R  radius  of  rotor 

r  distance  of  blade  element  from  axis  of  rotation 


x  abscissa,  along  axis  OX  of  the  rigid  blade,  of  an  element  of 

the  elastically  deflected  blade,  with  the  flapping  axis  as 
origin  (fig.  l) 

|  abscissa  along  axis  OX 

y  ordinate  of  a  point  of  the  elastically  deflected  blade  along 

_ an  axis  OY  normal  to  axis  OX _ 

iThe  "mode  0”  is  that  which  corresponds,  for  a  hinged  blade,  to  a 
vibration  without  bending,  that  is,  to  flapping  3.  The  mode  1  for  such 
a  blade  is  then  that  which  corresponds  to  a  vibration  of  one  node. 

^This  damping  term  has,  obviously,  no  significance  when  there  is 
separation  of  flow,  because  beyond  angle  of  separation  the  normal  lift 
coefficient  is  practically  independent  of  the  angle  of  attack. 
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distance  of  the  flapping  hinge  from  the  axis  of  rotation  of 
the  rotor 

flapping  angle  of  the  rigid  blade 

mass  of  the  blade  per  unit  length  at  a  point  under  consideration 
blade  chord 

Young's  modulus  of  the  blade  material 

moment  of  inertia  of  a  normal  section  of  the  blade 

angular  velocity  of  the  blade 

time 


azimuth  angle  of  the  blade  at  time  t 

phase  difference  of  the  azimuth  angle  due  to  damping 

forward  speed  of  helicopter 

resultant  velocity  of  the  air  on  a  section  of  the  rigid  blade 
normal  component  of  velocity  v 
tangential  component  of  velocity  v 
angle  of  attack  of  a  blade  section 

angle  of  the  resultant  velocity  v  with  the  normal  plane 


angle  of  velocity  v  with  the  resultant  velocity  on  the 


deflected  blade 


ry  cos  cp 


angle  of  the  normal  plane  with  the  forward  velocity 


ratio  of  advance 


air  density 

lift- curve  slope  of  the  profile 
acceleration  of  gravity 


k 
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static  moment  with  respect  to  the  axis  of  rotation  of  the  part 


of  the  blade  located  beyond  the  abscissa  x 


(f-; 


*  (a  +  |)d|) 


vi,o> 


coefficient  of  damping  due  to  deflection  of  the  blade.  (See 
text. ) 

natural  deflection  function  of  order  i 

natural  frequency  of  the  order  i  of  the  nonrotating  blade 

natural  frequency  of  the  order  i  of  the  blade  caused  by 
rotation  with  angular  velocity  to 


auxiliary  functions,  dependent  on  time  only 


n  order  of  a  harmonic  in  Fourier  series 

na,  nb  subscripts  of  terms  in  cosine  and  sine  of  a  Fourier  series 


THE  EQUATION  OF  DEFLECTION  OF  THE  BLADE 


Let  us  consider  the  forces  which  act  at  S  on  a  blade  element  of 
span  dx  in  the  plane  YOX.  (See  fig.  1.)  These  forces  are: 

(1)  The  elementary  lift3  dF  =  |  l  ^  ^ WT 2  -  VNVT^dx 

(2)  A  corrective  term  of  the  damping  of  the  lift  force  due  to  the 

flexural  elastic  deformation  of  the  blade.  This  term  is  of 
the  form 

-Ky  dx 

(3)  The  weight  dp  =  m’g  dx,  the  components  of  which  along  OX  and 

OY  are,  respectively 

m’g  sin  p  dx  (negligible) 
m'g  cos  p  dx  ~  m'g  dx 


is  elementary  lift  already  includes,  according  to  the  definition 
of  Vjj,  a  term  due  to  the  aerodynamic  damping  of  the  flapping  rigid  blade. 
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(4)  The  centrifugal  force  of  rotation  m'aPr  dx,  the  components  of 

which  along  OX  and  OY  are,  respectively 

m’ou£(a  +  x  cos  (3 )dx  ~  m'cu^r  dx 

and 

m'o3^(a  +  x  sin  (3)dx  o*  m’oo^rP  dx 

(5)  The  force  of  inertia  of  flapping:  -m'xp  dx 

(6)  The  force  of  inertia  of  deflection:  -m’y  dx 
Strictly  speaking,  the  following  forces  should  also  he  included: 

(7)  The  centrifugal  force  of  flapping:  m'xP2  dx 

(8)  The  Coriolis  force  due  to  the  simultaneous  action  of  blade 

deflection  and  flapping  motion: 

-2(3 yin'  dx 

however,  we  shall  disregard  them  in  relation  to  the  centrifugal  force  of 
rotation. 

For  the  blade  element  d|  to  be  in  equilibrium,  it  is  necessary  to 
add  to  these  forces  the  unknown  actions  of  the  adjacent  elements  on  the 
element  under  consideration. 

Calculation  of  Ky 


VT  our  +  V  cos  a  sin 
cos  cp  cos  <p 


£$  =  i-  C-°S  9  (See  fig.  2.) 
v 


y  cos 


our  +  V  cos  a  sin 
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It  becomes 


K  =  2.  i  ^12.  (usr  +  V  cos  a  sin  \|r) 
2  di 


Since  the  calculation  can  be  made  only  when  the  coefficient  of  the 
damping  term  due  to  the  deflection  is  independent  of  time,  the  periodic 
part  in  V  sin  \|/  will  be  disregarded.  This  reduces  the  problem  to  the 
corresponding  mean  speed  of  the  air  with  respect  to  the  blade  element. 


Hence 


n  dC7 

K  =  %  l  - £.  (1ST 

2  di 


The  differential  equation  of  the  blade  then  reads 


dxc 


f El  ^ 

—  f  m' (a  +  |)d| 

\  dx 

2 J  dx 

dx  J 

X 

;  +  m'y 


+  Ky 


HF  .  ..  .2 

—  -  m'g  -  m’xp  -  m'o)  rp 
dx 


(1) 


The  dots  signify  the  derivatives  with  respect  to  time. 


This  equation  of  the  partial  derivatives  must  be  completed  by  four 
limiting  conditions  which  define  the  integration  constants,  namely. 


y  =  0  for  x  =  0 
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METHOD  OF  SOLUTION 


The  foregoing  deflection  equation  is  composed  of  a  first  member  with 
terms  dependent  on  the  deflection,  and  of  a  second  member  which  is  inde¬ 
pendent  of  deflection.  This  second  member  comprises  the  aerodynamic 
forces,  the  weight  and  the  forces  of  inertia  of  rotation,  and  flapping 
acting  on  the  rigid  blade.  These  forces  are  easily  computed  for  each  point 
of  the  blade  and  for  each  one  of  its  azimuth  positions  after  the  equation 
of  flapping  p  of  the  rigid  blade  has  been  solved.  We  shall  waive  their 
calculation  and  identify  the  second  member  of  the  preceding  equation  by 
the  function  F'^(x,t).  We  further  put 


-R 


m 


(a  +  i)d! 


s 


so  that  the  preceding  equation  reads 


f?  (EI  §)  ‘  *  E  (s  &)  +  m'»  +  (2) 

This  is  an  equation  of  partial  derivatives  with  second  member, 
representing  the  forced  vibrations  of  the  blade  with  damping. 

To  resolve  it,  we  introduce  the  natural  functions  of  the  deflections 
of  the  blade,  that  is,  the  vibrations  obtained  by  solving  the  foregoing 
equation  of  the  partial  derivatives  above  without  second  member. 

Consider  first,  for  simplicity,  a  state  of  forced  vibrations  with¬ 
out  damping  arising  from  the  deflections  of  the  blade.  (In  other  words, 
the  term  Ky  is  disregarded. ) 

When  we  consider  the  moment,  Mrig,  of  the  forces  exerted  on  the 
assumedly  rigid  blade,  this  moment  is,  at  a  point  of  the  blade,  a  func¬ 
tion  of  time  only  and  can  therefore  be  developed  in  series  of  periodic 
functions  of  ^  =  cut 

Mrig  =  Mq  +  Ma  cos  t  +  Mb  sin  t  +  M2a  cos  +  Mpb  sin 


+  Mga  cos  3^  +  sin  3^  +  •  •  • 


8 


NACA  TM  1312 


This  moment  is  none  other  than  the  moment  of  the  distributed  outside 
forces  F*£  appearing  in  the  second  member  of  the  preceding  equation. 

Let  us  resolve  these  forces  at  each  instant  in  series  of  distribu¬ 
tions  such  that  they  each  impart  to  the  blade  a  deflection  taking  the 
form  of  a  natural  mode  of  deflection  of  the  corresponding  order. 

It  can  be  shown  that  an  arbitrary  deflection  of  the  blade  may  always 
be  resolved  in  series  of  natural  functions  by  reason  of  the  relations  of 
orthogonality  existing  between  the  natural  functions  of  continuous  beams, 
whatever  their  limiting  conditions.  (Physically,  it  means  that  the  dif¬ 
ferent  natural  vibrations  act  independent  of  each  other  without  mutual 
interactions. ) 

The  proposed  resolution  has  the  form 


F'd  =  +  +  82m*Tl2  +  •  *  •  (3) 

where  g^  is  a  function  of  time  only  and  is  the  natural  function 

of  the  deflection  of  the  blade  of  the  order  i. 

The  deflection  of  the  order  i  of  the  blade  affects  then,  at  each 
instant,  the  form  of  the  function  q^,  that  is,  it  will  be  given  by 


*i  =  hiTli 

where  h^  is  a  function  of  the  time  only.  The  total  deflection  is,  by 
virtue  of  the  relations  of  orthogonality  invoked  above,  obtained  by 
superposition  of  the  various  natural  deflections  of  the  blade  vibrating 
at  the  corresponding  natural  frequencies,  with  amplitudes  and  phase  dif¬ 
ferences  defined  by  the  function  h^,  which,  itself,  is  obtained  by 
putting  the  expression  y.^  =  h^q^  in  the  equation  of  the  partial 
derivatives. 

Hence,  for  a  frequency  of  the  order  i 

hl  d?  (EI  S1) ' hia?  (s  sr) +  =  gi'°’’'i 


(M 


NACA  'EM  1312 


9 


Now  it  is  known  that  an  equation  of  the  form 


» 


+  m*y  =  0 


representing  the  natural  vibrations  of  a  freely  vibrating  blade  and 
involved  in  rotation  with  an  angular  velocity  ou  permits  an  infinite 
number  of  solutions  of  the  form  q  sin  vt,  satisfying 


and  the  limiting  conditions.  Every  solution  is  the  natural  func¬ 

tion  of  order  i  corresponding  to  the  natural  frequency  ^ 

Therefore,  when  q^  is  a  natural  function,  it  simultaneously 
satisfies  the  equation  (4)  and  the  equation 


v2i,oin’Tli 


where  is  the  corresponding  natural  frequency  of  the  blade  actuated 

with  a  speed  of  rotation  cu. 

Hence,  after  simplification 

hiv2i,u)  +  *1  =  gi 

Since  the  functions  hj_  and  gj_  are  periodic  with  respect  to  \Jr, 
they  can  be  developed  in  harmonic  series 

hi  =  hi jO  +  hij a  cos  \|r  +  hi,b  sin  ^  +  hi, 2a  cos  2\Jr 

+  hi^b  sin  2^  +  hi,3a  cos  +  hi^b  sin  3^  •  •  • 

«i  =  Si,0  +  Si, a  cos  *  +  Si,b  sin  *  +  Si, 2a  cos  2^ 


+  Si, 2b  sin  2^  +  gi, 3a  cos  3^  +  §i,3b  sin  3^  .  . 
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The  differential  equation 


therefore  resolves  itself  hy  identification  of  the  coefficients,  and 
limiting  i.t  to  the  third  harmonic,  it  reduces  to  the  system 


5 

i,<x> 

Si,0 

hi,a( 

,2  - 
1,00 

■  J)  = 

Si, a 

*2  . 
1,03 

-2)  = 

Si,  b 

hi,2a( 

'v2i,- 

-  4<a2) 

=  Si,  2a 

hi,2b< 

(v2i,0) 

-  4a)2) 

=  Si,  2b 

hi,3a' 

3 

•v 

•H 

OJ 

-  9“>2) 

gi,3a 

hi,3bi 

(y2i,o) 

-  9a)2) 

=  Si,  3b 

Each  function  gi  is  well  defined. 

To  determine  it,  simply  multiply  the  two  members  of  equation  (3) 
by  tj-l  and  integrate  over  the  blade. 


(5) 


Owing  to  the  conditions  of  orthogonality 
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it  leaves 


R 


R 


F'd^i  **  =  gi  m,Tli 
0  J0 


»«  2dx 


hence 


R 


'0 


F* d^i  ** 


Si  = 


(6) 


-R 


>0 


T  2, 

m1 1|^  dx 


The  bending  moment  exerted  on  the  elastically  deformable  blade  is 
computed  next. 

The  elastic  deformation  due  to  a  single  harmonic  of  the  outside 
forces,  such  as 


1=00 


i=l 


is,  as  shown  previously. 


y>~  (gi,nam,TU  cos  n*  +  gi,nbm,TU  sin 


1  =  00  1  =  00 

i=l  i=l  1 


On  replacing  the  terms  h^  by  their  values  obtained  from  equation  (5), 
the  corresponding  bending  moment  reads 


1=00 


M 


elast. 


•  -  k  -  X 


i=l 


g 


i,n  [d  d2q,  cos 

El  1  J 


X  -  ifiaP  dx2  1  sin  n\)r 
1;U)  L 


(7) 
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On  making  the  exact  calculation  of  the  natural  functions  t]^,  it  is 
found  that  they  vary  very  little  with  a>  and  that  a  natural  function 
of  the  order  i,  for  00  =  0,  can  be  compared  with  the  same  natural 
function  for  the  normal  speed  o>. 


''i. 


(!) 


This  simplification  is  not  exact,  but  the  resulting  error  is  small 
(less  than  3  percent  for  the  first  natural  function  in  the  case  of  the 
blade  cited  in  the  example  hereinafter) ,  being  of  the  same  order  as  ' 
those  committed  in  the  distribution  of  the  masses  and  the  flexural 
stiffness  of  the  blade. 

In  this  case,  the  natural  function  satisfies  both 


Hence 


El 


dx  d! 


On  replacing  El  -  by  its  value  in  (7),  the  coefficient  of  the 

dx^ 

harmonic  n  of  the  elastic  moment  which,  in  fact,  bends  the  flexible 
blade,  is  given  by  the  expression 
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The  value  gj_  is  computed  by  (6)  for  the  values  of  \j |r  for  which  F*  ^ 

is  given.  The  development  of  in  Fourier  series  defines  the 

coefficient 


Si,n 


The  natural  functions  q^  and  the  integrals  JJ' m'q^  dx  d|  can  be 

computed  by  a  classical  method  such  as  the  iteration  method  (Stodola), 
or  the  Galerkin  method. 

There  is  no  occasion  to  be  preoccupied  with  normalizing  the  natural 
functions.  For  the  amplitude  of  q^,  any  convenient  scale  is  suitable; 

the  effect  of  the  scale  disappears  later  in  the  product 


«i  fj  m'ni  a*  <U 


The  calculation  for  <n  =  0  is  made  while  remembering  the  previous 
statement  that  the  deflection  is  practically  unmodified  by  rotation. 

If  the  iterative  method  is  used  for  computing  the  natural  functions 
of  the  deflection,  the  natural  frequencies  of  the  blade  not  rotating  and 
of  the  blade  rotating  at  angular  velocity  o>  can  be  computed  by  applying 
Rayleigh's  method  to  the  obtained  natural  deflection.  This  method  affords 
rigorous  solutions,  converges  rapidly,  and  avoids  the  solution  of  n  equa¬ 
tions  with  n  unknowns  to  which  the  Galerkin  method  leads. 

Note:  The  bending  moments  could  also  be  computed  by  direct 
application  of  equation  (7).  This  method  is  predicated  on  the 
exact  knowledge  of  d^^/dx^  which  prohibits  the  use  of  the  approxi¬ 
mation  =  q^  q  because  a  slight  error  in  a  function  can  cause 

a  substantial  error  in  its  second  derivative.  The  calculations  of 
the  natural  deflections  are  quite  complicated.  On  the  other  hand, 
the  function  y  being  defined  by  dots,  it  is  not  possible  to  derive 
it  directly  to  obtain  d^y/dx^. 
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SIMPLIFICATION  OF  THE  CALCULATIONS 


Numerous  calculations  made  on  blades  of  various  helicopters  with 
different  plan  forms  and  distributions  of  masses  and  different  amounts 
of  rigidity  have  shown  that,  in  certain  cases  -  blades  of  moderately 
conical  shape,  little  twist,  and  lightly  loaded  at  the  tip  -  the  natural 
functions,  other  than  the  first,  exert  little  influence  on  the  maximum 
bending  moments  exerted  on  the  blade,  and  consequently  on  the  maximum 
alternating  fatigue  stresses  to  which  the  blade  is  subjected.  This  is 
because  the  distributions  of  the  outside  forces  (curves  F*d  =  f(r)  for 

the  various  i)  for  such  blades  represent  the  behavior  of  the  first 
natural  distribution  (curve  m'q  =  f (r) )  in  a  satisfactory  manner  and 
also  because  the  natural  functions  of  higher  orders  present  all  the 
"loops"  and  "nodes"  in  continuously  increasing  number,  matching  poorly 
the  greater  deflection  which  the  maximum  moments  produce. 

It  is  only  in  cases  of  small  deflections  that  the  "parasite  remainder" 
of  the  higher  frequencies,  arising  from  the  fact  that  the  distributions 
of  the  outside  forces  never  have  a  curve  exactly  identical  with  the  first 
natural  distribution,  can  play  a  significant  part. 

When  the  blades  have  pronounced  camber  and  twist  and  are  lightly 
loaded  at  the  tips,  the  forces  of  inertia  can  become  more  important  than 
the  aerodynamic  forces  at  the  blade  tip  in  the  entire  sector  of  the  swept 
disk  corresponding  to  the  advancing  blade.  The  distributions  of  the  out¬ 
side  forces  F'd  =  f(r)  can  assume,  therefore,  the  curves  approached  by 
the  second  natural  distribution  (curve  m'r^)  for  an  entire  series  of 
and,  in  that  case,  the  second  and  sometimes  the  third  natural  functions 
must  be  taken  into  consideration  in  the  calculation  of  the  maximum 
moments  exerted  at  the  blade. 

Examination  of  the  curves  F'd  =  f(r)  permits  one  to  determine, 
with  a  little  practice,  when  resolution  of  the  outside  forces  can  be 
limited  to  the  first  natural  distribution  and  the  deflection  of  the 
blade  to  that  corresponding  to  the  first  natural  function. 

In  the  latter  case,  the  calculations  are  considerably  simplified. 

The  distribution  of  the  outside  forces  is  reduced  to 


F'd  =  g1m*Tj1 
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The  tending  moment  on  the  elastic  blade  becomes 


M 


elast.  n 


1,0) 


-  n2cs2 


gl,n 


,R 


m'T)^  dx  d| 


(9) 


On  comparing  this  bending  moment  with  the  moment  of  the  outside 
force  distribution  F’^  exerted  on  the  assumedly  infinitely  rigid  blade, 
we  get 


M 


rig 


If  F'd  ^  di  =  gl 


m*T]^  dx  d| 


(10) 


It  is  seen  that  the  bending  moment  on  the  elastic  blade  is  obtained 
by  multiplying  the  coefficients  of  the  harmonics  of  the  same  order  of 
moment  on  the  rigid  blade  by  a  factor 


An  - 


-1*0. 


v  ^ 

1,0) 


(11) 


which  is  constant  over  the  blade  for  a  given  harmonic. 

It  is  no  longer  necessary  to  calculate  the  natural  functions 
and  the  natural  frequencies  V-j,  q  and  ^  They  can  be  readily  and 

closely  approximated  by  the  Rayleigh  method  applied  to  a  curve  repre¬ 
senting  approximately  the  deflection  of  the  blade  while  still  satisfying 
the  limiting  conditions  rigorously. 

As  regards  the  various  harmonics  to  be  kept  for  the  moments  on  the 
rigid  blade,  it  seems  that  no  advantages  are  gained  by  going  beyond  the 
third,  which  is  already  relatively  small. 

The  bending  moment,  at  a  point  on  the  abscissa  x,  on  the  rigid 
blade  is  then  given  by  the  expression 


Mrig  =  MO  +  Ma  cos  \Jr  +  M-jj  sin  +  M£a  cos  + 


Mgb  sin  2\| 1  +  M^a  cos  3i  +  sin  3t 


1 6 
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The  effective  moment  bending  the  actual  elastic  blade  at  the  same 
point  is  obtained  by  the  simple  relation 


v2i,o 

Melast.  =  — -  M0  + 


v2i,o 


V^i  -  (X? 

1,03 


^Ma  cos  \|r  +  M-jj  sin  + 


- — : — -  (M?a  COs  +  M?b  sin  2k)  + 

1,03 


~1~" -  (M3a  cos  ^  +  M3b  sin  3t) 


1,03 


-  9032 


CALCULATION  INCLUDING  AERODYNAMIC  DAMPING  DUE  TO 
DEFLECTION  -  GENERAL  METHOD 


In  this  case,  the  term  Ky  disregarded  in  equation  (2)  must  be 
included. 

It  has  been  shown  that 


OIT 


The  term 


is  homogeneous  to  a  distributed  mass. 


Putting 
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gives 


K  =  <D(x)m,a) 

dCz 

where  \  is  a  dimensionless  coefficient.  Since,  as  a  rule,  £■  l  — —  r 
(x)  2  di 

is  not  proportional  to  the  distributed  mass  m' ,  the  coefficient  ®(x) 
is  a  function  of  x. 


Having  already  permitted  one  approximation  for  K  by  including 
only  the  corresponding  mean  speed  on  a  blade  element,  another  simplifica¬ 
tion  is  effected  by  bringing  only  the  mean  value  of  along  the 

span  of  the  blade  into  the  first  member  of  the  equation  of  motion. 


d>  = 


R  -  a 


p  i£^dx 

2  m«  di 


Accordingly,  equation  (2)  reads 


IL/ei  a?  3-  (a  $iL)  +  m’y  +  tm'tty  =  F*  (12) 

dx2\  d x2/  dxV  W  d(x,r) 


With  the  symbols  already  employed,  we  get 


hiv  i,(s  +  hi  +  =  «i 

The  foregoing  system  of  equations  (5)  is  therefore  replaced  by 
another  of  the  form 


hi,na(v2i,oo  "  n2^)+  hi,nb^ 
hi,nb(v2i,cs  -  n2^-  hi,nana)2<f 


®i,na 

Si,nb 


(13) 
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Consequently 

hi,na  = 


(Ha) 


The  calculation  is  then  carried  out  in  the  same  way  as  without 
damping;  the  deflection  due  to  a  single  harmonic  n  of  the  outside 
forces  being  always 


i=oo 


Replacing  hj_  by  its  value  gives 


Hence,  as  before,  the  term  of  the  harmonic  n  of  the  moment  which,  in 
fact,  bends  the  flexible  blade 


«i,nb  sin  ^i,n)cos  n*  +  (&i,na  sin  *i,n  +  «i,nb  cos  *i,n)sin 
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or  also 


If,  as  previously,  the  distribution  of  the  outside  forces  is  limited 
to  the  first  distribution  and  the  deflection  of  the  blade  to  that  corre¬ 
sponding  to  the  first  natural  function,  the  following  rule  results:  The 
effective  bending  moments  acting  on  an  elastic  blade,  with  damping,  are 
obtained  from  the  bending  moments  acting  on  the  assumedly  rigid  blade, 
each  harmonic  of  which  is  modified  as  follows: 

(l)  Its  amplitude  is  multiplied  by  a  coefficient  equal  to 


(2)  There  is  a  forward  phase  difference  of  \|/n,  so  that 


2 

n$o> 


tan  = 
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PRACTICAL  PROCEDURE  OF  CALCULATION  IN  THE  CASE  OF  THE  GENERAL  METHOD 
(The  Double  Resolution  of  the  Outside  Forces  Being  Limited  to  the 
Third  Natural  Distribution  and  to  the  Third  Harmonic  in  >10 


The  first  step  is  to  determine  the  natural  functions  of  the 
blade  and  =  f*( r) ,  as  well  as  the  corresponding  natural 

frequencies  Vj_  q ,  Vg  q ,  and  0  nonrotating  blade  and  the 

V1  ay  v?  ay  an<^  cd  blade  rotating  at  angular  velocity  u>. 

y  y  y 


Next,  it  is  necessary  to  evaluate  the  quantities 


and  then  plot  the  curves 


JJ  m*Tj1  dx  d|. 


rrR 

/  /  m'qg  dx  d|,  and 


dx  dl 


against  r. 

Then  determine  the  outside  forces  on  the  rigid  blade  F’^  =  f(r) 
for  different  >jf  (eight  positions  spaced  k^°  apart  must  be  explored). 

Next  evaluate  for  each  position 


Thence  one  obtains 


g3  = 


L  r’a''3ax 


m*q. 


"dx 


for  each  position 
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Develop  g^,  g^,  g^  in  Fourier  series  (by  the  Runge  method,  for 

example)  and  stop  with  the  third  harmonic  of 

whence 


;1,0 

gl,a 

gl,b  * 

'  ’  gl,  3"b 

;2,0 

S2,a 

g2,b 

S2,3b 

3,0 

g3,a 

g3,b  * 

*  *  S3,3b 

Calculate  the  bending  moments  on  the  blade  for  different  stations 
x  =  K.R  by  the  formula 


i=3 


i=l 


with 
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v2i,a>  "  n2a^ 

To  determine  quadrants  of  \|r,  <p:  For  ^  n,  -  >  0  and 


for  cpi?n,  g;L,na/coS  <pi#n  >  0. 


cob  *i#n 


PRACTICAL  EXAMPLE 


Hereinafter  follovs  an  application  of  this  method  to  the  calculation 
of  the  hending  moments  exerted  on  the  blade  of  the  helicopter  N.C.  2001  in 
forward  flight,  with  p  =  0.43. 

The  mechanical  characteristics  of  the  blade,  which  has  a  radius  of 
R  =  6.85m,  are  given  in  figure  3*  Figures  4  and  5  give  the  natural 
functions  q  and  the  blade  distributions  m'q.  Figure  6  gives  the 
natural  frequencies  of  the  nonrotating  blade  and  of  the  blade  rotating 
at  angular  velocity  cn  for  the  different  modes  of  vibration.  Figure  7 
gives  the  distribution  of  the  aerodynamic  forces,  and  figure  8  the  dis¬ 
tribution  of  the  total  outside  forces  F’^  on  the  rigid  blade  for 
several  azimuth  positions  of  the  latter. 

According  to  figure  8,  the  distributions  of  the  outside  forces  for 
the  positions  cp  =  45°,  90°,  135°,  and  l80°  have  clearly  the  shape  of 
the  second  natural  distribution  m’T^j  which  explains,  as  will  be  seen 
in  figures  9  to  12,  the  importance  of  the  bending  moments  computed  with 
the  second  natural  function  included. 

Figure  13  shows  the  bending  moment  distribution  for  several  radii 
plotted  again  \lr ,  and  figure  l4  the  enveloping  curve  of  the  maximum 
bending  moments  exerted  on  the  blade.  Figure  14  also  shows  the  maximum 
bending  moment  curve  acting  on  the  blade  of  the  S.E.  3000  helicopter 
(R  =  6m,  p  =  0.4l),  whose  mechanical  characteristics  are  shown  also  in 
figure  3.  It  is  seen  that  for  this  blade,  which  presents  an  average 
camber  and  is  lightly  loaded  at  the  tip,  the  error  made  in  the  maximum 
bending  moment  by  limiting  the  distributions  of  the  outside  forces  to 
the  first  distribution  does  not  exceed  6  percent,  which  justifies  the 
simplification  of  the  calculation  indicated  previously. 

For  the  blade  of  the  NC.2001,  the  maximum  bending  moment  is  severely 
subjected  to  the  influence  of  the  second  distribution,  but,  by  way  of 
compensation,  it  is  practically  clear  of  that  of  the  third  distribution 
(except  at  the  tip). 

However,  it  should  be  noted  that  the  latter  assumes  a  significant 
part  in  the  evaluation  of  the  negative  maximum  bending  moment,  expressed 


2b 
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by  Mmi- n  (see  figs.  10  and  13) ,  and  consequently  in  the  appraisal  of 
the  maximum  alternating  fatigue,  defined  by 


f  =  + 
max 


^max  “  ^min 
2W 


where  W  =  ^  is  the  resistant  modulus  of  the  particular  section. 


Translated  by  J.  Vanier 
National  Advisory  Committee 
for  Aeronautics 
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Figure  3.-  Plan  forms  and  distribution  of  mass  and  rigidity  of  the  S.E.3000 

and  N.C.2001  helicopter  rotor  blades. 


Figure  4.-  N.C.2001  blade  -  natural  elastic  deflection  curves  of  the  blade. 
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Figure  8.-  N.C.2001  blade  -  total  forces  F'^  distributed  over  the  rigid  blade 


I 


M(mkg; 


Figure  11.-  N.C.2001  blade  -  blade  bending  moments  for  i|r  =  180°. 
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Figure  12.-  N.C.2001  blade  -  blade  bending  moments  for  *  =  270°. 
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First  natural  mode  only 
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modes  included 
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